Effect of Non-Magnetic Impurities (Zn, Li) 
in a Hole Doped Spin-Fermion Model for Cuprates 
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The effect of adding non-magnetic impurities (NMI), such as Zn or Li, to high-Tc cuprates is 
studied applying Monte Carlo techniques to a spin-fermion model. It is observed that adding Li is 
qualitatively similar to doping with equal percentages of Sr and Zn. The mobile holes (MH) are 
trapped by the NMI and the system remains insulating and commensurate with antiferromagnetic 
(AF) correlations. This behavior persists in the region %NMI > %MH. On the other hand, when 
%NMI < %MH magnetic and charge incommensurabilities are observed. The vertical or horizontal 
hole-rich stripes, present when % NMI=0 upon hole doping, are pinned by the NMI and tend to 
become diagonal, surrounding finite AF domains. The %MH-%NMI plane is investigated. Good 
agreement with experimental results is found in the small portion of this diagram where experimental 
data are available. Predictions about the expected behavior in the remaining regions are made. 



PACS numbers: 74.62.Dh, 74.20.Mn, 71.10.Fd 
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Neutron scattering studies of the high-Tc cuprates 
have established that spin incommensurability (IC) ap- 
pears in these compounds upon hole doping, fil There 
is mounting evidence indicating that this IC is due to 
the formation of charge stripes. |2j One way of under- 
standing the relationship between magnetic and charge 
fluctuations, as well as their impact in other properties 
of these materials, is by introducing NMI in the Cu-0 
planes. Some of these experimental results are puzzling. 
It has been observed that replacing Cu 2+ by Zn 2+ de- 
pletes antiferromagnetism without causing spin IC even 
at 25% doping. It is estimated that an impurity con- 
tent as large as 30-40% will be required to destroy AF 
order in this case. p]|| While the addition of MH (Sr in 
La 2 - x Sr x Cu04 (LSCO)) to lightly Zn-doped materials 
produces spin IC similar to the one observed in LSCO, 
H doping with Li, which introduces equal numbers of lo- 
calized impurities and mobile holes does not induce spin 
IC at 10% Li doping of La 2 Cu0 4 (LCO). § 

The goal of this paper is to gain theoretical under- 
standing on the effect of NMI doping on cuprates. Nu- 
merical studies of Hubbard and t-J models, the Hamil- 
tonians traditionally used to study the physics of the 
cuprates, are very difficult at low temperatures. A sim- 
pler spin-fermion model has been used successfully to un- 
derstand magnetic and charge properties in these mate- 
rials, showing that spin IC and charge stripe formation 
are related. || The spin-fermion model is constructed 
as an interacting system of electrons and spins, crudely 
mimicking phenomenologically the coexistence of charge 
and spin degrees of freedom in the cuprates. B-lljj]. Its 



Hamiltonian is given by 



H 



(ij)" 



h.c. 



J^s; -Si + J'^Si-Sj, 

i <ij> 

w 

where c\ a creates an electron at site i = (i x , i y ) with spin 
projection a, s\=}2 a p c\ a a a pC\p is the spin of the mo- 
bile electron, the Pauli matrices are denoted by a, S; is 
the localized spin at site i, (ij) denotes nearest-neighbor 
(NN) lattice sites, t is the NN-hopping amplitude for the 
electrons, J > is an AF coupling between the spins of 
the mobile and localized degrees of freedom, and J' > 
is a direct AF coupling between the localized spins. The 
density (n)=l — x of itinerant electrons is controlled by a 
chemical potential /i. Hereafter t = 1 will be used as the 
unit of energy. ^From previous phenomenological analy- 
sis the coupling J is expected to be larger than t, while the 
Heisenberg coupling J' is expected to be smaller. jR||ll],[| 
The value of J will be here fixed to 2 and the coupling 
J' = 0.5, as in Ref. ||. To simplify the numerical cal- 
culations, avoiding the sign problem, localized spins are 
assumed to be classical (with |Si|=l). This approxima- 
tion is not as drastic as it appears and it was discussed 
in detail in Ref. ||. The model will be studied using 
a Monte Carlo (MC) method, details of which can be 
found in Ref. ]12|]. Periodic boundary conditions (PBC) 
are mainly used but it has been checked that the results 
presented here are independent of the boundary condi- 
tions. 

Doping with Zn will be simulated by randomly remov- 
ing localized spins. This is achieved by turning J = 
at the impurity site and J' = along the four links that 
connect the impurity to neighboring sites. Li doping has 
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to be treated differently. Since it is in an oxidation state 
Li + (while Cu and Zn are 2+), the ion effectively acts 
as a negative charge center. To simulate this effect the 
hopping t that connects the impurity to neighboring sites 
should be reduced. In addition, the Li ion is closed shell 
and thus, non-magnetic. This feature is simulated as in 
the case of Zn, by randomly removing localized spin de- 
grees of freedom. An additional effect of Li doping is 
the introduction of one MH, thus the chemical potential 
has to be tuned such that the percentage of mobile holes 
equals the percentage of localized impurities. 

To study the magnetic properties of the system the 
spin-spin correlation functions among the classical spins 
defined as u>(r) = ^ ■ Si+ r ) (where N is the number 

of sites) was measured. Its Fourier transform provides the 
spin structure factor S(q). The momentum q 7 takes the 
values 27rn/L 7 , with n running from to L 7 — 1, and L 7 
denoting the number of sites along the 7=x or y direction. 
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FIG. 1. (a) Structure factor 5(q) for the localized spins 
versus momentum, for 3—2, J'=0.05 and T=0.01 on an 
8x8 lattice and several Zn concentrations at (n)=l along 
(0,0) — (0, 71") — (n, tt) — (0,0); (b) spin and charge distribu- 
tion for a typical MC snapshot for 15% Zn and the same 
parameters as in (a). The arrow lengths are proportional to 
the projection of the localized spin on the (X-Z) plane, and 
the radius of the circles is proportional to the local electronic 
density n(i), according to the scale shown. PBC are used, 
(c) Density of states for different percentages of Zn doping at 
(n) = 1 with the same parameters as in (a). The dark areas 
indicate states with very little dispersion 

In Fig.l-a the structure factor for different concentra- 
tions of localized NMI, like Zn, is shown for an 8 x 8 lat- 
tice at temperature T = O.Oli and at (n) = l since doping 
with Zn does not modify the density of itinerant elec- 
trons. It can be seen that the peak remains at (tt, it) for 
all the studied NMI concentrations. The only effect of 



adding NMI is to reduce the intensity of the structure 
factor since spins are being removed. It is clear that ran- 
dom impurities do not destroy the AF order as can be 
observed from a snapshot at 15 % doping for the spins on 
the plane X-Z (Fig.l-b). The snapshot also shows a uni- 
form charge distribution. This is in complete agreement 
with experimental data for Zn doping in LCO. || Our 
simulations can produce dynamical information directly 
in real-frequency without the need of carrying out (un- 
controlled) analytic extrapolations from the imaginary 
axis. This allows us to study the behavior of the den- 
sity of states (DOS). In Fig.l-c the DOS as a function of 
%NMI doping is presented. The addition of localized im- 
purities creates states inside the gap which are indicated 
by a black filling in the figure. Analyzing the spectral 
functions A(q, to), it was observed that these states have 
very small dispersion which added to the absence of spec- 
tral weight at the chemical potential indicates that the 
system remains an insulator, in agreement with experi- 
ments, n 
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FIG. 2. (a) Schematic magnetic and charge phase diagram 
in the plane %Sr — %Zn. The thick lines and the dark areas 
indicate the regions that have been studied experimentally; 
(b) Spin and charge distribution for a MC snapshot on an 8x8 
cluster for 3=2, J' =0.05, T=0.01, at (n)» 0.85 and % Zn=0. 
The spins are shown on the X-Y plane; (c) same as (b) but 
with % Zn=5; (d) DOS for different cases with %Sr>%Zn. 
The notation is as in Fig.l. PBC are used 

The next step in our investigation is to study magnetic 
and charge properties varying both the mobile hole (Sr) 
and localized NMI (Zn) densities. In this case only a 
small fraction of the %Sr-%Zn plane has been explored 
experimentally. Studies were performed for < %Zn < 
25_and < %Sr < 3 @, < %Zn < 75 and %Sr = 15 

1, < %Zn < 12 and 15% Sr 
I Sr Bil, < %Zn < 40 and 



12% Zn and 14% Sr { 
< %Zn < 4 and 10' 
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0% Sr |J. Numerous experiments have also been per- 
formed on LSCO (without Zn). For < %Sr < 25 neu- 
tron scattering experiments are presented in Ref. [ jl6| . 
Zn-doped YBa2Cu307_s (YBCO) has also been studied 
with < %Zn < 4 and the O content varying between 6.6 
and 7, which roughly corresponds to a percentage of MH 



ranging between 9-15%. |Ti 



15| . These regions are shown 



in Fig.2-a. In the same figure the magnetic and charge 
properties of the model Eq.(l) in the plane %Sr — %Zn 
obtained from the present study are presented. In the 
absence of NMI, the addition of MH induces stripe for- 
mation, as well as spin IC. The stripes are either hor- 
izontal or vertical || and a snapshot for 15% MH and 
no localized impurities is shown in Fig.2-b. As it can be 
seen in the figure, the electronic density is depressed near 
the stripe indicating that perpendicular fluctuations are 
present, i.e., there is short range movement of holes in 
the direction perpendicular to the stripe. 

When localized NMI are added they tend to trap holes. 
This reduces the effective hopping towards the impurity's 
nearest neighbor sites and diagonal stripes anchored by 
the impurities tend to develop as boundaries of AF do- 
mains. This effect is similar to the formation of diago- 
nal stripes due to pinning experimentaly observed. JlS| ] 
An example of this behavior is shown in Fig.2-c where 3 
NMI (i.e. 5% Zn) have been added to the system with 
(n) ~ 0.85. A 7T— shift also occurs between the AF do- 
mains separated by the diagonal stripes. For some distri- 
butions of Zn impurities both diagonal and non-diagonal 
stripes, pinned by the impurities, coexist. 
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FIG. 3. (a) Snapshot for %Zn=%Sr=15 in the plane X-Z; 
(b) S(q) for different points with %Zn > %Sr; (c) Density of 
states for the parameters shown in (b). The parameters and 
notation are as in Fig.l 

It has also been observed that while a localized im- 
purity does not disturb the AF background, when the 
number of mobile holes is larger than the number of local- 
ized impurities a n~ shift occurs in the x and y directions 



across the impurity which eventually will give rise to spin 
IC (see Fig.2-c) @|o). According to this result, spin IC 
should be detected in the region %Sr > %Zn, more specif- 
ically, when the percentage of Sr is larger than that of Zn 
by approximately 5% according to experimental results 
for LSCO. This is in agreement with our observations if 
it is taken into account that due to the finite size of our 
system, spin IC is first detected when the peak in the 
structure factor is at (it, 37r/4) and symmetrical points. 
Thus, it is here predicted that spin IC should be exper- 
imentally observed in samples with %Sr > %Zn. Notice 
that, due to the fact that the stripe structure is modi- 
fied by NMI doping this may not be the same kind of 
IC observed in LSCO. In fact, it is possible that besides 
the magnetic peaks at (w — S, 7r) and symmetrical points, 
less intense peaks may appear at (n — S, 7r — 5). More 
generally, as the percentage of NMI increases, IC peaks 
maybe observed in all directions in momentum space. 
The DOS indicates that in this situation the system has 
metallic characteristics although the chemical potential 
is in a pseudogap as shown in Fig.2-d. This feature is also 
in agreement with experiments, pi] ] A study of the spec- 
tral function shows that the states indicated with black, 
have very little dispersion. From the figure it can be seen 
that some mobile holes are trapped in these states and 
become almost localized but the extra mobile holes go 
into the dispersive states that form the pseudogap. 
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FIG. 4. (a) Density of states for different percentages of Li 
doping on an 8x8 cluster for J=2, J'=0.05, T=0.01;(b) Struc- 
ture factor for the localized spins as a function of the momen- 
tum. The parameters are the same as in (a); (c) snapshot for 
%Li=15 (same parameters as in (a) and same notation as in 
Fig.l.) 

When %MH < %NMI, the spinless centers tend to 
localize the mobile holes as it is shown in Fig.3-a for 
%Zn = %Sr = 15. Though the hopping towards the im- 
purity sites is t = 1, the electrons prefer to stay away from 
that site, since the coupling J is zero and the observed 
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electronic density in them is just (ni) « 0.03, which in- 
dicates that the holes are very localized at the impurity 
sites. Notice that the hole localization, which is observed 
experimentally appears naturally in our calcula- 

tions, while it has to be assumed in other theoretical 
approaches. The system is magnetically commensu- 
rate and the structure factor peaks at (n, tt) as shown in 
Fig.3-b for different values of %Zn > %Sr. In these cir- 
cumstances the system is an insulator as it can be seen 
in Fig.3-c where the DOS is displayed. When %Zn=%Sr 
the chemical potential lies in a gap separating dispersive 
from non-dispersive states. When %Zn > %Sr the chem- 
ical potential is located at the non-dispersive states. In 
both cases the system is an insulator. As mentioned be- 
fore, very few experiments where performed in this region 
of the plane. However, our results are in agreement with 
the available data since insulating behavior was found 
for %Sr = 15 and %Zn > 15 Ref. g| and %Sr = 3 and 
%Zn > 3 Ref. (§. 

Note that the behavior observed along the line 
%Sr = %Zn is in qualitative agreement with the exper- 
imental results available for Li doped LCO. It is well 
known that the system is an insulator up to the maxi- 
mum possible percentage of Li (50%) p3| , p4j and at 10% 
doping spin IC is not observed Q which was very unex- 
pected. Thus, these results indicate that qualitatively, Li 
doping should be equivalent to equal percentages of Zn 
and Sr doping. However, in order to represent Li doping 
more accurately, it has to be considered that Li is less 
positively charged than Zn and Cu and thus it should 
act as an attractor of mobile holes. To mimic this be- 
havior the hopping towards the impurity sites should be 
reduced. In our previous results it has been observed 
that t=l tends to greatly localize the mobile holes since 
the electronic density at the impurity sites is only 0.03. 
Thus, any smaller value of t will not introduce important 
qualitative changes but, for completness, and since it has 
been used in other approaches, J22| the results for t=0 
near NMI are also presented. In this case the MH get 
trapped by the impurity and become totally localized. 
This can be seen in Fig.4-a where the density of states 
is shown for different percentages of Li dopping. Local- 
ized states (marked in black), which do not disperse with 
momentum, proportional to the number of added holes 
appear in the density of states, and the chemical potential 
lies inside a gap indicating that the system is insulator as 
experimentally found. p5| , p4| . Incommensurability is not 
observed as it can be seen in Fig.4-b where the structure 
factor is displayed. A snapshot is presented in Fig.4-c. 

Summarizing, the effect of adding non-magnetic im- 
purities to a model of interacting fermions and classical 
spins has been investigated using MC techniques. When 
the percentage of mobile holes is larger than the NMI, 
charge and spin IC is observed. The NMI act as pin- 
ning centers for the horizontal and vertical stripes that 
are observed in the absence of vacancies. The pinned 
stripes tend to become diagonal and surround AF do- 
mains. Thus, the observed incommensuration may not 



be identical to the one observed in LSCO. When the per- 
centage of MH is equal to or smaller than the NMI, the 
mobile holes are trapped by the impurities and the sys- 
tem remains insulator and magnetically commensurate. 
Thus, based on our study, it is predicted that incommen- 
suration should not be experimentally observed in this 
situation, and that samples with equal percentages of Sr 
and Zn should behave similarly to the Li doped ones. 
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